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ABSTRACT 

Extensive spectral observations of -q Carinae over the last cycle, and particularly 
around the 2003.5 low excitation event, have been obtained. The variability of both 
narrow and broad lines, when combined with data taken from two earlier cycles, reveal 
a common and well defined period. We have combined the cycle lengths derived from 
the many lines in the optical spectrum with those from broad-band X-rays, optical 
and near-infrared observations, and obtained a period length of Pprcs = 2022.7 ± 1.3 
d. 

Spectroscopic data collected during the last 60 years yield an average period of 
Pavg = 2020 ± 4 d, consistent with the present day period. The period cannot have 
changed by more than AP/P — 0.0007 since 1948. This confirms the previous claims 
of a true, stable periodicity, and gives strong support to the binary scenario. We have 
used the disappearance of the narrow component of He i 6678 to define the epoch of 
the Cycle 11 minimum, Tq =JD 2,452,819.8. The next event is predicted to occur on 
2009 January 11 (±2 days). The dates for the start of the minimum in other spectral 
features and broad-bands is very close to this date, and have well determined time 
delays from the He i epoch. 

Key words: stars: general - stars: individual: eta Carinae - stars: binary. 



1 INTRODUCTION 

T] Carinae is one of the most luminous stars in the Milky Way 
and contains many mysteries. It has been attracting atten- 
tion since the 1820s, when it suffered large brightness fluc- 
tuations, culminating with the giant eruption that ejected 
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the Homunculus in 1843. The star faded to naked eye in- 
visibility, and after the discovery of the supernovae in the 
XlXth century it was classified as a slow supernova. How- 
ever, around 1940, it started to brighten again, indicating 
that the star was only hidden by dust, not destroyed. 

The spectrum is rich in emission lines of low excita- 
tion species: Hi, F en, [Feiil, [Ni ii], Tin, etc. (|Thackeravl 
Il953l ): after 1944 (|Gaviolal Il953 ) high excitation forbid- 
den lines of [Neni], [Arm ], [Siii], and [Feni] c an also be 
readily identified (see also Damineli et al.l [19981 . and refer- 
ences therein). Today we know that the narrow lines (for- 
bidden and permitted) are emitted in the Weigelt blobs 
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l|Weigelt fc Ebersberged 198d), at ~ 0.3 arcsec from the 
central star ( Davidson et alJ Il995h and the broad emis- 
sion hnes are formed i n the wind of the c entral object 
jHillier and Allenlll992l : [Davidson et alJ 119951 ). The combi- 
nation of high and low excitation lines in the same object, 
however, was paradoxical. 

A key to understanding this interesting object was 
found recently through the study of the variability of 
the high excitation lines. The high excitation forbidden 
lines disappeared in 1948, and again in 19 65, 1981 , 1987 
and 1992. These 'sp ectroscopic event s' dCaviola 1 19531 : 
iRodgers fc Searldl 19671 : [Thackcrav 1967: Zanella et al.lll98j ') 
or 'low excitation events' (|Damineli et al.l Il998t ) were be- 
lieved to be part of S Doradus cycles, commonly seen in 
other LBV stars similar to eta Car. This interpretation 
seemed to be supported by the He I A10830 line which went 
to minimum ( |Damin cli I ggj) when the near-in frared light 
curve went to maximum jWhitelock et al]|l994l '). The max- 
ima in the near-infrared light curves were not truly peri- 
odic and the length of the quasi-period was different for dif- 
ferent pass-bands. However, the spectrosco pic events were 
demonstrated to be periodic (|Damineli et al...200Q ). in con- 
trast to the incoherent character of the S Dor oscillations. 
iDamineli. Conti and Lopesl(| 19971 ) and others proposed a bi- 
nary model with a highly eccentric orbit, a hotter secondary 
component and a strong wind- wind collision (WWC). Bina- 
rity is interesting as it potentially allows the direct measure- 
ment of the mass of the stars, their most fundamental pa- 
rameter. The binary scenario has provided a framework for 
understanding the star and provided guidelines for fruitful 
observations, although so me prefer a model in which there 
are periodic shell ejections (jMartin et al.|[2006l ). In Fig. [T] we 
present examples of high and low excitation state spectra of 
Tj Carinae. 

The observation of an event in 1997.95, as was pre- 
dicted, brought more confidence to the t rue periodic nature 
of the variation (jDamineU et al.ll2000l ). iFeast et al] (120011 ) 
used archival spectra to identify three previously unreported 
events, in 1953, 1959 and 1970, which also fit the 5.5-yr 
period. Moreover, those authors discovered that the dips 
on top of the broad quasi-periodic near-infrared maxima 
were truly periodic and correlated with the behavior of the 
high excitation lines. An extensive X-ray monitoring cam- 
paign was started in 1996 with the RXTE satellite and re- 
vealed deep minima in 1997.95 and 2003.49 which coincide d 
with the minima seen at other wavelengths (|Corcoranll2005[ ). 
X-ray observations inside and outside the minimum per- 
formed with Chandra and XMM furnished details on the 
column density (A^h), temperatur e and chemical compos i- 
tion of the colliding wind shock (|Hamaguchi et al] [20o3) . 
Ivan Genderen et al.l (|2006h showed that the optical light 
cur ve displays periodic dips like those in the near-infrared 
and lLaiiis etail (|2003h reported a very detailed light curve 
in the B, V, R and I bands for th e 2003.49 event. The event s 
were recorded also at radio-cm (IDuncan and White! I2OO3I ) 
and radio-mm (| Abraham et aLll2005l ), but no specific value 
to the period length was reported for those wavelengths. 

Many other features vary periodically in intensity and 
radial velocity, like the broad emission and P Cygni ab- 
sorption components, and can also be used to derive the 
period length. One o f them is Hell A4686 discovered by 
ISteiner and Daminelil (|2004l ). which raises and drops just 



before minimum faster than any other eature over the en- 
tire spectrum. Although faint (EW< 2 A) it was frequently 
monitored with high signal/noise along the last event. Un- 
fortunately, it was observed only occasionally in the previous 
events, precluding its use to measure the period. This spec- 
tral line deserves better monitoring in future events, not only 
to improve the accuracy of the derived period, but also be- 
cause it is the highest excitation feature observed at optical 
wavelengths, and its origin remains a mystery. 

To facilitate di scussion we label the eve nts by num- 
bers as described bv lGroh and Daminelil (|2005l ) : number one 
(#1) is assigned to the event observed in 1948 by Gaviola, so 
that the event of 2003.49 is #11. We define cycle as the time 
interval between the starting of two consecutive minima, so 
that cycle #9 started at the 1992.42 minimum and finished 
when cycle #10 was starting in 1997.95. Because of observa- 
tional reasons, that will become clear later in this paper, the 
starting point of a cycle is defined by the disappearance of 
the He I A6678 narrow line component. With this definition, 
phases along the cycle are defined in a unique way for all 
measured quantities. 

The paper is organized as follows. We present in section 
(21 the observations; in 331 the definition of the phase of the 
minimum; in §3) the determination of the period length; in 
^ the stability of the period; in ^ the relation between 
the sharp peaks during the giant eruption and periastron 
passages; and in 33 the discussion and conclusions. 



2 OBSERVATIONS AND MEASUREMENTS 

The majority of the ground-based observations presented in 
this paper came from a monitoring campaign started in 1989 
at the Coude focus of the 1.6-m telescope of Pico dos Dias 
Observatory (OPD-LNA/Brazil). The observational setup 
at OPD was kept essentially unchanged through the cam- 
paign: a dispersion grating with 600 1/mm, entrance slit 
width '^1.3 arcsec, exposure time ~ 5 s in Ha increasing to 
~ 15 min at 3500 and 10800 A. Spectra were extracted along 
~ 2 arcsec in the spatial direction and no measurable differ- 
ences in line intensities were seen when changing the extrac- 
tion size by a factor of 2. Three different CCDs have been 
used, with resolving powers R = 25 km s~^ (0.25 A pixel"^) 
at Ha in 2003 and R = 50 km s"^ (0.39 A pixel"^) in the 
preceding years. On some occasions, a 1024 x 1024 Hawaii 
detector was used to observe the He I A10830 line, deliver- 
ing a spectral resolution R = 40 km s~^ (0.65 A pixel"^). 
On other occasions, spectra of this line were taken at R = 
15 km s~^ with a thinned CCD. After correcting for fringes 
and degrading the spectral resolution, these spectra were 
almost identical to those collected with the infrared array 
at the same date. Observations in 1992 and 1997/8 were 
done with a thick CCD that was almost free of fringes, but 
had a low sensitivity in the blue, which explains the poor 
coverage of important lines in that spectral range. For wave- 
lengths longer than 6500 A, telluric absorptions and fringes 
(in thinned CCDs) were removed by using templates con- 
structed from spectra of bright early type stars {9 Carinae, 
C Ophiuci or Puppis) observed immediately after or before 
77 Car. 

For the 2003 event, we also used spectra taken with the 
spectrographs REOSC (R = 25 km s"^) and EBASIM (R = 
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Figure 1. Sample of the rj Car spectra at high and low excitation states, (a) high resolution red spectrum showing the disappearance of 
narrow emission lines and strengthening of P Cygni absorption components during the minimum, (b) low resolution blue spectra showing 
the high state in 1995 and two spectra during the minimum taken in 1997 December 31 and 1970 May 17. 
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Figure 2. Definition of phase - (a) He I A6678 line profiles observed along 17 days in 2003, showing the disappearance of the narrow 
component, (b) The fading phase, showing our method to derive the phase (JD 2,452,819.8) of He I A6678 and A7065 narrow line 
components. 



7 km s ^) attached to the 2.15-m CASLEO telescope (Ar- 
gentina), and spectra taken at CTIO with the 4.0-m Echelle 
Spetrograph (R = 8 km s"^) and at Magellan with MIKE 
Spectrograph (R = 12 km s"^). For the 1997/8 event, we 
used spectra collected at La Silla/ESO with CAT-CES (R 
= 12 km s"'^). For the 1992 event, we also used spectra 
collected with the FLASH/HEROS spectrograph attached 
to the 50-cm telescope (ESO/Chile) with a fiber diameter 
~ 5 arcsec and spectral resolution R = 12 km s~^. On sev- 
eral occasions we used FEROS spectrograph attached to the 
1.52-m telescope at La Silla to cover the entire optical win- 
dow at resolution R = 12 km s~^. 

Before measuring the spectral features, we degraded the 
spectra to a dispersion of 0.39 A pixel"^. This step was 
not really necessary but it helped facilitate the adoption of 



the same limits between the narrow and broad components, 
and positioning of the stellar continuum, when measuring 
the spectra. Since we adopted the observations collected at 
LNA Observatory as a reference, we added data from other 
sources only in the case where they merged smoothly to the 
line intensity curve. This criterion was fulfilled by almost 
all ground-based observations, confirming our expectation 
that slit widths in the range 1 — 3 arcsec width would give 
the same results independent of the position angle (P.A.) of 
the slit. This happens because the main emitting region is 
smaller than 1 arcsec and has a huge contrast to the sur- 
rounding Homunculus nebula and also because the seeing 
fwhm is larger than 1 arcsec, smearing out the emitting re- 
gion. In a forthcoming paper (on the long term behavior of 
the spectral lines) we will present the complete list of obser- 
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vations from the entire campaign and a table with individ- 
ual measurements. Fig. [1^ displays spectra representative 
of the high and low excitation states, showing the disap- 
pearance of the high excitation lines and enhancement of 
P Cygni absorption profiles during the minimum. Fig. [ijj 
shows spectra in the blue for the high excitation state of 
1995 and for the low excit ation state of 1997 and 1970 (see 
also iDamineli et al.l (|l998t ) for the full spectral range 3850- 
11000 A). 

Spectra from the Space Telescope Imaging Spectrome- 
ter (STIS) on the Hubble Space Telescope are available for 
the 2003.49 and 1997.95 events, though for consistency we 
do not include them here since the slit width is much nar- 
rower than used in the ground-based observations, sampling 
only a part of the inner circumstellar nebulosity. These data 
are of course important for disentangling stellar from cir- 
cu mstellar varia t ions, a n d have been more ful l y described 
inlNielsen et al.1 (l2007al'). iNielsen et al.1 (|2007bh . IGuU et al.1 
(|2006l ). and [Davidson et al. I (|2005^ . Since the w nd of the 



primary star is resolved by the STIS slit and the slit's posi- 
tion angle varied in different visits, care must be taken when 
comparing line profiles from different epochs. This applies to 
the lower excitation transitions, formed far from the central 
source(s) that may be subject to spatial asymmetries. 

AU the data processing and measurements were done in 
the standard way using IRAF packages. Narrow lines were 
modeled by Gaussian fitting and deblended from the broad 
components. Since they are seated on top of broad line pro- 
files, which are themselves variable, we referred their equiv- 
alent widths (EW) to the local stellar continuum, in order 
that these measurements correspond to line flux normalized 
to the local stellar continuum, instead of classical equivalent 
width. As in the case of EW, this kind of measurement is 
translated into line flux when multiplied by the stellar con- 
tinuum flux. Because of this, we use the simple designation 
of equivalent width in place of normalized line flux. Broad 
line emission profiles were separated from the narrow com- 
ponents, when they existed, and their equivalent widths and 
baricenters (for radial velocities) were measured by direct in- 
tegration along the line profile. Radial velocities are in the 
heliocentric reference system. 

It is difficult to attribute errors to single measurements, 
as the main source is systematic, not statistical. The spectra 
were well exposed, in order that photon noise is very low, 
except in the violet region. The major source of error is 
linked to the stellar continuum, because of line blendings 
and changes in relative intensity of line/continuum, as the 
seeing changes and smears out the central source of emission 
lines. The random errors can be judged by the smoothness of 
the curves in line intensity and the plots show that they are 
small, in general comparable to the size of the symbols in the 
figures. We minimized the errors by over-plotting the spectra 
and pointing the cursor always in the same postilion. We 
must warn, however, that this procedure does not eliminate 
the systematic errors. 



3 DEFINING PHASE FOR THE 
SPECTROSCOPIC EVENT 

A simple method to measure the periodicity of the events 
is through the disappearance of spectroscopic features like 



the high ionization lines or the narrow components of He I 
(Fig. [2^). In practice this is difficult because of the follow- 
ing: a) the time sampling has been too coarse to pick up the 
exact time when the feature disappears; b) spectroscopic 
features reach minimum at different times; c) minima are 
usually reached asymptotically for many important spectral 
features, often taking up to a week to disappear completely; 
d) when the line EWs are less than ~ 100 mA, they are dif- 
ficult to measure, unless the spectra have very high signal- 
to- noise ratio (S/N). In addition to producing a large un- 
certainty in the epoch of the minimum, faint features may 
not be directly connected to the emitting region, but can 
be light echoes that fall inside the slit aperture. Moreover, 
in some cases a very faint blended line remains in emission 
through the minimum, as in the case of [Arm] A7135. 

In order to minimize these problems we restricted our 
analysis to the phase of steep decline, which lasts for ~ 2 
weeks, starting ~ 3 weeks before complete disappearance. 
We performed a linear fit, and extrapolated it to zero in- 
tensity to determine the time of minimum (Fig. [2Jd) . This 
procedure is much more robust than other techniques, since 
it does not require a dense time sampling along the mini- 
mum. It is relatively insensitive to the S/N of the spectrum, 
and is easily reproducible by other observers. The epoch of 
minimum, i.e. phase (the starting point of the deepest part 
of the minimum) for the He I A6678 narrow line component 
derived by this method is To = JD 2,452,819.8 (2003 June 
29 or 2003.491). Since He I A6678 is strategic for the spectro- 
scopic event (it has a long observational history, shows clear, 
easily measured variability and lies in a spectral range with 
good CCD efficiency) we chose it for our definition of phase 
0. 

There are two situations for which it is useful to find 
signatures that indicate the time of phase 0: when examin- 
ing non-calibrated historical spectra or when trying to track 
the evolution of an event during a monitoring campaign. As 
the high excitation lines are much more variable than the 
lower excitation lines, and because the spectrum has plenty 
of lines, it is relatively easy to find line pairs that interchange 
peak intensity ratio with time. A high excitation line, as the 
minimum approaches, decreases until its peak is equal in 
strength to that of some nearby low excitation line (in gen- 
eral Fell or [Fell]), and we record the date when this occurs. 
The faster the high excitation line varies the more accurate 
is the determination of the time of change in the line ratio. 
This happens for dates close to phase zero, when the vari- 
ability is high, but we were able to find good line pairs up 
to three months before phase and almost two years after. 

We display in Table [1] the time in days for the inversion 
in peak intensity ratios, relative to phase 0. Negative val- 
ues represent dates before phase and positive values dates 
later than phase 0. Entries in column 2 are for the fading 
phase and in column 3 for the recovering phase, except for 
He I A10830, which displays the two ratio inversions in the 
fading phase. Times are shorter in column 2 than in column 
3 due to the fact that the fading phase is fast and the recov- 
ering phase is slow. From an examination of data for the last 
three cycles, we found that times in column 2 are accurate 
to ~ 15 per cent and in column 3 to ~ 25 per cent. 

In Table[T]we have also listed the line He I A10830. It has 
a double peak, like in classical Be stars. The V ('violet') and 
R ('red') peaks are variable, both in intensity and in their 
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Figure 3. Period length: (a) from the EW He I A6678 narrow line component along the latest 2 events; (b) from radial velocity of the 
He I A6678 broad component in the last 3 events; (c) from the /C-band by using the method of folding the 2 latest events and minimizing 
the residuals. 



Table 1. Time delays in days, relative to phase 0, when the 
intensity of line peaks change ratio. 



Line ratio 


change to < 1 


change to > 1 


He I 10830 R/V* peaks 


-105 


-4 


He: 4471/ [Fell] 4475 


-46 


+550 


[Sm] 6312/Fen 6317 


-15 


+442 


[Fein] 4658/ [Fell] 4640 


-17 


+358 


He I 7065/ [Fen] 7171 


-9 


+148 


He I 6678/ [Ni 11] 6666 


-9 


+145 


[Felll] 4658/ [Fell] 4475 


-8 


+250 


[Feiii] 4658/ [Fell] 4665 


-5 


+168 


He: 5876/Nai 5890 





+18 


[Nil] 5754/ [Fell] 5746 


+1 


+79 


Fell 8490/Feii 8499 


+7 


+40 



* 'red' (R) and 'violet' (V) peaks intensity ratio 



relative strength. For almost the entire 5.5-yr cycle, R > V. 
As the minimum approaches, the R peak starts decreasing 
faster than the other, in such a manner that 105 days before 
phase they reach V = R, changing to R < V subsequently. 
The R < V state lasts for almost three months when the 
rate of fading of the R peak slows down and the V peak 
starts falling fast. Just four days before phase 0, the peaks 
again reach V=R and return to R > V. 



4 THE PERIOD LENGTH 

There are a number of ways to measure the period length; 
the best one for spectroscopic data is based on the He I nar- 
row line components. The equivalent width of this feature 
remained relatively constant at ~1500 mA for most of the 
cycle #10. About three weeks before phase 0, it began to 
change fast, declining by ~ 25 mA day^^. We used this 
segment of the line intensity curve to measure the period, 
applying a scheme of epoch foldi ng and minimiza tion of dif- 
ferences similar to that used by ICorcora 3 (|2005h . Since we 
sampled better the fading phase to the minimum, it was 
sufficient to shift this piece of the line intensity curve from 
the event #10 until it matched that of event #11 (Fig. [3^) 



Table 2. Period length in days from different spectral regions 



Period ± error spectral feature or pass-band 



2026 


± 2 


He I 6678 narrow component Eq. W. 


2024 


± 2 


X-rays 


2023 


± 1 


J-band 


2023 


± 1 


_ff-band 


2023 


± 1 


ii'-band 


2023 


± 2 


L-band 


2022 


± 2 


Fe II 6455 P Cyg abs. radial velocity 


2022 


± 1 


Si 11 6347 Equivalent Width 


2022 


± 1 


He I 6678 broad radial velocity 


2022 


± 1 


He I 10830 Equivalent Width 


2021 


5 ± 2 


y-band 


2021 


± 2 


Fe II 6455 P Cyg abs. Eq. W. 


2022 


7 ± 1.3 


average ± standard deviation 



to derive the period. We repeated the same procedure with 
the event #9, getting the best fit for P = 2026 days with 
an uncertainty of 2 days. A careful examination of Fig. [3^, 
however, indicates that the slope of the fading phase was 
steeper during event #10 than in event #11. This is due to 
secular changing in the line intensity and this is the main 
source of errors in the period determination by this method. 

The broad component of He I A6678 also changes quickly 
before the minimum. Its radial velocity decreases slowly in- 
between the events, but three weeks before phase it re- 
verses the trend, and starts to increase. About 4 weeks after 
phase 0, the radial velocity increases at a rate of ~ 5km s~^ 
per day. The steep rise in the radial velocity curve is useful 
to determine the period length, in the same way as we have 
done for intensity of the narrow component. This procedure 
might be more robust than using equivalent widths, since 
radial velocities are much less affected by the secular vari- 
ations in intensity. By combining RVs from the last three 
events (#9, #10 and #11) we derive P = 2022 ± 1 day (Fig. 

Other spectral features were observed, whenever possi- 
ble, and some of them also were useful for measuring the 
period length. The total equivalent width of the He I A10830 
line recorded in the last three events gives P = 2022 ±1 days. 
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The equivalent width of the P Cygni absorption component 
of the Fe ll A 6455 A line (also using the last three events) 
gives P = 2021 ± 2 days and the radial velocity curve of the 
same line (events #9 and #11) results in P = 2022 ±2 days. 
The equivalent width of Si ll A 6347 A P Cygni component 
(events #9 and #11) resuhs in P = 2022 ± 1 days. 

These periods ar e in good agreem ent with P = 2024±2 d 
derived from X-rays (jCorcoranlpOOsI'l. with P = 2023 ± 3 d 
from the K-h&nd light curve JWhitelock et al.ll2004h. and 
with P = 2021.5 ± 0.5 d from V-band (|van Genderen et all 
I2OO6I ). We re-derived the periods obtained by those au- 
thors, using their published data and applying the tech- 
nique of folding the pre- minimum branch of two events, 
#10 and #11 in the cas e oflWhitelock et all (|2004ll and #7 
and #11 in the case of Ivan Genderen et al.1 (|2006i 'l. In the 
case of F-band photometr y, we combine d the light curve o f 
Ivan Genderen et al.1 lj2006l l with that of lLaiiis etHI (j2003l l 
in order to get a better definition of the descending branch 
of the 2003.5 minimum. The derived peri od was the same 
as published bv lvan Genderen et al.l (j2006h . but with an un- 
certainty of 2 d instead of 0.5 d. In the case of near-infrared, 
IWhitelock et al.l (|2004[ ) used the lower point in the J^-band 
minimum. Our procedure of minimization of residuals ap- 
plied to the JHKL—ha.nd photometry (Fig. [S];) gave the 
same period, but with a tighter constraint. It is encouraging 
to see that the period is robustly defined, independent of 
the particular choices made by different authors during the 
measurements (Table [2}. 

Since there is no reason to suppose that the period 
length would depend on the particular technique used, we 
combined all these individual periods to get a mean value 
to the period. We call it the present day period (Ppros) to 
differentiate from that determined from historical observa- 
tions. Since the systematic errors may be more important 
than statistical errors, we report the uncertainty in the pe- 
riod as a simple standard deviation. 

Ppros = 2022.7 ± 1.3 d or Pp^os = 5.538 ± 0.004 yr. 

Regarding the times of phase 0, there is no reason to 
expect that different features give the same epoch, since they 
are produced in a variety of regions - in the stellar winds, 
in the WWC and in the circumstellar material. Since we 
are dealing mostly with spectroscopic lines we define, for 
reasons discussed earlier, phase of the spectroscopic events 
from the disappearance of the He I narrow line intensity. This 
yields the following ephemeris: 

JD(phase 0) = 2, 452, 819.8 ± 0.5 + (2022.7 ± 1.3d)E. 

The uncertainty is only 0.07 per cent of the period 
length, which enables us to accurately predict the time of 
phase of the next spectroscopic event: JD 2,454,842.5 ± 2 
(2009 January 9-13). 



5 STABILITY OF THE PERIOD 

An important question is the long term stability of the pe- 
riod, since the companion stars are losing mass at high rates 
and tidally interact during the periastron passages. In ad- 
dition, the primary star could be a fast rotator - as in- 
dicated by its dense polar wind (jvan Boeckel et al] l2005l : 



Table 3. Period derived from epochs of predicted and observed 
minima. 



Cycle 


Pred. 


Predicted 


Observed 


P 




JD (**) 


Date 


Date 


days 


#1 


32,592.8 


1948 Feb 11 


1948 Apr 19 


>2015.9 


#2 


34,615.5 


1953 Oct 26 


1953 Jun 28* 


<2029.3 


#2 






1953 Dec 30 


> 2009 


#3 


36,638.2 


1959 Mar 11 


1959 May 14 


>2014.4 


#4 


38,660.9 


1964 Sep 23 


1965 Feb 15 


>2001.9 


#5 


40,683.6 


1970 Apr 07 


1970 May 17 


>2016.0 


#7 


44,729.0 


1981 May 04 


1981 May 21 


>2018.6 


#8 


46,751.7 


1986 Nov 17 


1987 Jan 15 


>2003.1 


#9 


48,774.40 


1992 June 01 


1992 May 31 




#10 


50,797.10 


1997 Dec 15 


1997 Dec 12 




* Intermediate excitation, ** JD-f- 


2,400,000 





IWeigelt et"al]|2007l : [s"mith et al.ll2003l ) - and the specific an- 
gular momentum may be changing continuously. The cou- 
pling of rotational and orbital angular momentum may lead 
to an increase in the orbital period. A period derived from 
data encompassing many cycles may be hiding such varia- 
tions. The only way to tackle this question is by compar- 
ing the present day period (measured from the last few 
events) with the average period (Pavg), derived from the 
events recorded in the past 60 years. 

We looked at the events reported bv lFeast et al.l (|200lh . 
In case the spectrum was in low excitation state, we assumed 
that it had just reached the minimum and so was at phase 
0. As a matter of fact, phase must have occurred already 
at some time before the observation, in order that a period 
based on that record would underestimate the period. The 
average period, or more precisely a lower limit to it, is ob- 
tained from the time interval between that date and the last 
recorded minimum (2003.49) divided by the number of cy- 
cles. The derived average periods are presented in column 
5 of Table [S] We excluded the last three events, since they 
were used to derive the present day period. The dates of 
observed m i nima in Table [3] (column 4) were taken from 
iFeast et al.1 (|200ll '). except for the observation on 1970 May 
17. This spectrum (taken at CTIO) was recorded by Virpi 
Niemela and indicates that the phase occurre d at least 8 
days b efore Thackeray's observation reported bv lFeast et al.l 
(|200ll ) . The CTIO logbook reports that spectra were taken 
by Barry Lasker a day before phase of the 1970 event 
(April 6) , which could lead to a very tight constraint on the 
period, but unfortunately we were not able to locate that 
spectral plate. 

The observation made on 1948 April 19 gives Pavg > 
2015.9 d and that of 1981 May 21 gives Pavg > 2018.6 d. 
We can constrain the period length also from the other side. 
A maximum period may be derived when a particular ob- 
servation was made before phase 0. This is the case for the 
observation made on 1953 June 28, when the star was ap- 
proaching the minimum, but was still in an intermediate 
phase, which gives Pavg < 2029.3 d. The average period is 
thus constrained to: 



The stability of the period can be obtained from the 
difference between the present day period and the average 
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Figure 4. Relative intensities of [Nil] A5754 and [Fen] A5746 - (a) Line profile variations. Labels at right of [Nil] peak indicate days 
from phase of the 2003.49 event, (b) Ratio between line peaks, indicating that [Nll]/[Fell]< 1.0 starts at phase and lasts for 78 days 
during the minimum. 



period, taking into account that Pavg refers to half of the 
cycles involved. The spectrum taken in 1953 June 28 indi- 
cates that the period cannot have decreased by more than 
1.4 d cycle"^ and that of 1948 April 19 implies that it cannot 
have increased by more than 1.4 d cycle"^. 

We have another way to constrain the average period 
using quantita tive information of the first event in 1948. 
iGaviolal (119531 ) reported that [N ll] A5754 was fainter than 
[Fell] A5746 , which places the date of the observation in 
a particular range inside the low excitation event. Exami- 
nation of recent events indicates that before phase 0, [N ll] 
is much stronger than the neighboring [Fe ll] line. The [N ll] 
line decreases quickly, in contrast to [Fe ll] which undergoes 
small and slow changes. Both features reach equal intensity 
0.7 days after phase 0, as can be seen in Fig,|4^, where varia- 
tions during event #11 are displayed. The [N ll] line remains 
fainter than [Fe ll] for a subsequent 78 days. This can be seen 
in Fig. |4}3, which combines measurements made in the last 
three events. The ratio of these two lines is not sensitive to 
the slit width or to the spectral resolution, as long as they 
are kept < 4 arcsec, or R > 2000, respectively. The fact that 
the 1948 observation was done < 78 days later than phase 
0, combined with the epoch of the 2003.49 minimum results 
in: 

Pavg = 2020 ± 4 d. 

This period is compatible with that derived in the 
present day data. Taking into account that the average was 
taken between 10 cycles, it could not have changed by more 
than 1.5 d cycle"^, in close agreement with the 1.4 d cycle"^ 
derived before. The constraint to the period change is: 

-0.0007 < AP/P < +0.0007. 

This c onfirm s previous claim s of st rict pe riodicity by 
iDaminelil (119961 ): iDamineh et~aLl (|2000l ) and iFeast et all 



(|200ll ) implying that the low excitation events are only un- 
derstandable in the binary scenario. It is also consistent 
with the period change expected because of mass-loss from 
the primary star. Simple considerations show t hat P/P = 
aMJM where a is a constant of order unity l|Khaliullinl 
1 19741 ). Ignoring changes in eccentricity we find: 

P = o.lla ( -\ days/cycle. 

\^10-3Moyr-i )\ M J ' ' ^ 

This is luUy consistent with the observed upper limit. 



6 WERE THE PEAKS IN THE GIANT 

ERUPTION PRODUCED BY PERIASTRON 
PASSAGES? 

IDaminelil (jl996h pointed out that the three most pronounced 
peaks observed in 1827-1843 were in close coincidence with 
the predicte d times of phase 0, when using the period 
P = 2014 d. iFrewl (|2004h also noted that the 'Lesser Erup- 
tion' which began on 1887.5 was within a few months of 
phase 0. When using the new period, derived in this work 
(P = 2022.7 d), the 1827.087 peak is now at phase 0.15; that 
of 1837.967 is at phase 0.13: and that of 1843.3 at phase 0.05 
(|Frewll200j ). while the 1887.5 event corresponds to a phase 
of 0.96. The correlation between peaks and the start of the 
spectroscopic events worsens with the new ephemeris. How- 
ever, an exact correlation between the times of the peaks and 
phase is not really relevant, since the time sampling of the 
visual light curve was not very dense and the real maxima 
could have been missed by the observers. Moreover, if the 
mechanism that produced those peaks is the same as the one 
that produces the broad max;ima observed presently in the 
near-infrared light curve, the lack of coinci dence with phase 
would not be a surprise. As reported bv IWhitelock et al.l 
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the JHK- light curves present maxima around phase 
0, but not in exact coincidence. These near-infrared peaks 
are quasi-periodic and may remain in high state for up to 
3 years, depending on the wavelength. Although the near- 
infrared light-curve is anti-correlated with the radio flux, it 
still can be explained as free- free emission if optical depth 
effects are taken into account (|Whitelock et al]|2004l ). 

In the l/-band, the maxima associated with the perias- 
tron are inconspicuous, as compared to those in the near- 
infrared. Could they have been more pronounced during the 
giant eruptions of the nineteenth century? There is no rea- 
son to believe so, as Car was much brighter in the optical 
than it is now, diminishing the contrast between the quasi- 
periodic maxima and the underlying stellar light. However, 
this is unsafe terrain, as we do not know what mechanism 
produced the pronounced peaks during the giant eruption. 

Recent estimates of the Homunculus mass suggest that 
more tha n 10 Mp, of mater ial was ejected during the great 
eruption (|Smith et al.ll2003l ). Given the larger mass-loss that 
occurred in the giant eruption it is likely that the orbital pa- 
rameters changed substantially during the great eruption, 
and thus it is not surprising there is not a one to one corre- 
spondence between orbital phase and the pronounced peaks 
observed in 1827-1843, even if they were a result of a binary 
interaction near periastron. 



7 DISCUSSION AND CONCLUSIONS 

We presented a homogeneous set of spectra covering the 
events #9 (1992.42), #10 (1997.95) and #11 (2003.49). We 
derived the period by measuring intensities of narrow lines 
(Weigelt blobs) and broad emission lines (stellar wind) and 
radial velocity variations from broad line components. These 
data, and others collected from the literature, enable an ac- 
curate determination of the period Ppros = 2022. 7±1. 3 d. 
An average period encompassing the past eleven cycles was 
found to be Pavg = 2020± 4 d, compatible with the present 
day period. The period change is smaller than 1.5 d cycle"^ 
along the last half century. 

It is difficult to imagine any mechanism other than or- 
bital motion which could maintain such a high stability, 
while allowing individual features to show the distinct light 
curves that are observed. No luminous unstable star is ex- 
pected to follow such a precise clock. Even if it is conceivable 
that a shell ejection could be involved in the periodic events, 
it should be triggered by the periastron passage, when the 
secondary star approaches the primary to a few stellar radii. 
However, in an incoming paper we show that the event starts 
when the secondary star is still near apastron. 

We must expect that the period is drifting, since the 
stars are losing mass, they interact strongly as they get very 
close at periastron and the primary star appears to be a 
fast rotator. The period changes have been AP/P < 1/1000 
along the last 60 years, which is consistent with the observed 
mass loss rate. 

The situation is different for the eruptions of 1843 and 
1890, when a considerable amount of matter was removed 
from the primary star in brief episodes. However, the re- 
lation between the pronounced peaks observed during the 
great eruption and periastron passages continues to be un- 
known. This is because the coincidence is not perfect, and 



because the peaks might not be strictly periodic but could 
still be associated with periastron passages, as seen presently 
in the near-infrared light curves. 

From the disappearance of the He I A6678 narrow com- 
ponent we determined the epoch of the start of the mini- 
mum to beTo = JD 2,452,819.8. The procedure to define the 
minimum requires fitting and extrapolating the line inten- 
sity variation along the descending part of the line intensity 
curve, in the two weeks preceding the minimum intensity. 
Because of this definition, and since the starting time of the 
minimum is different from line to line, this definition is arbi- 
trary and has no physical meaning. However, it is robust and 
demands only a few observations along the ~3 weeks before 
the complete disappearance of the feature. Importantly, the 
time delay for all other features to reach the minimum is 
well known. 

The next minimum is predicted to start on 2009 Jan- 
uary 11 (±2 d). This will the best event since 1948 for 
ground-based observations, since its central core fits entirely 
in the good observing season. The next favorable event will 
not occur before 2020. In order to improve the results pre- 
sented in this work, daily observations should be made along 
a month starting on 2008 December 20. 
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